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Abstract 

An effective Hamiltonian for the Bose system in the mixture of ultracold atomic 
clouds of bosons and fermions is obtained by integrating out the Fermi degrees 
of freedom. An instability of the Bose system is found in the case of attractive 
interaction between components in good agreement with the experiment on the 
bosonic 87 Rb and fermionic 40 K mixture. 



The seminal paper by N.N. Bogoliubov "To the theory of superfluidity", published in 
1947, had a great influence on the development of the modern theory of quantum many- 
body systems. It is cited in all textbooks devoted to this topic. However, for a long time 
the real experimental systems which could be quantitatively described by the model of 
dilute Bose gas developed in this article were absent. Only about 50 years later, the Bose- 
Einstein condensation was achieved in the dilute ultracold gas of magnetically trapped 
alkali atoms. Since the first realization of Bose-Einstein condensation in ultracold atomic 
gas clouds [TJ E], studies in this direction have yielded unprecedented insight into the 
quantum statistical properties of matter. Besides the studies using the bosonic atoms, 
growing interest is focused on the cooling of fermionic atoms to a temperature regime 
where quantum effects dominate the properties of the gas [3]. This interest is mainly 
motivated by the quest for the Bardeen-Cooper-Schrieffer (BCS) transition in ultarcold 
atomic Fermi gasses 

Strong s-wave interactions that facilitate evaporative cooling of bosons are absent 
among spin-polarized fermions due to the exclusion Pauli principle. So the fermions are 
cooled to degeneracy through the mediation of fermions in another spin state [3J E] or 
via a buffer gas of bosons [S] (sympathetic cooling). The Bose gas, which can be cooled 
evaporatively, is used as a coolant, the fermionic system being in thermal equilibrium 
with the cold Bose gas through boson-fermion interaction in the region of overlapping of 
the systems. 

However, the physical properties of Bose-Fermi mixtures are interesting in their own 
rights and are the subject of intensive investigations including the analysis of ground 
state properties, stability, effective Fermi-Fermi interaction mediated by the bosons, and 
new quantum phases in an optical lattices. Several successful attempts to trap and cool 
mixtures of bosons and fermions were reported. Quantum degeneracy was first reached 
with mixtures of bosonic 7 Li and fermionic 6 Li atoms. Later, experiments to cool mixtures 
of 23 Na and 6 Li, as well as 87 Rb and 40 K jS], to ultralow temperatures succeeded. 
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In this article we study the instability and collapses of the trapped boson-fermion 
mixture due to the boson-fermion attractive interaction, using the effective Hamiltonian 
for the Bose system jHl Ej. We analyze quantitatively properties of the 87 Rb and 40 K 
mixture with an attractive interaction between bosons and fermions recently studied by 
Modugno and co-workers [3] . They found that as the number of bosons is increased there 
is an instability value Nb c at which a discontinues leakage of of the bosons and fermions 
occurs, and collapse of boson and fermion clouds is observed. Using the experimental 
parameters we estimated the instability boson number Nb c for the collapse transition as 
a function of the fermion number and temperature and found a good agreement with 
experimental results. 

First of all we briefly discuss the effective boson Hamiltonian [SJ U\ ■ Our starting point 
is the functional-integral representation of the grand-canonical partition function of the 
Bose- Fermi mixture. It has the form [HJ: 

Z= ( D[cj>*]D[(t>]D[r]Dty]exv{-\{S B {ct>\cl>)+ (1) 



h 

and consists of an integration over a complex field 0(r, r), which is periodic on the 
imaginary-time interval [0, h/3], and over the Grassmann field ip(r, r), which is antiperiodic 
on this interval. Therefore, 0(r, r) describes the Bose component of the mixture, whereas 
ip(r, r) corresponds to the Fermi component. The term describing the Bose gas has the 
form: 

" h<i r ( f d h 2 v 2 

S B {d>\d>) = I dr I dr\ <f (r, r)[h—-- + 
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+ V b (t)- ^)0(r,r) + ^|0(r,r)| 4 ). (2) 

Because the Pauli principle forbids s-wave scattering between fermionic atoms in the same 
hyperfine state, the Fermi-gas term can be written in the form: 

S F (r,i>) = / dr / dv{r(r,r) ih— + 



dr Imp 

+ Vp(r)-/^)V>(<r,r)}. (3) 

The term describing the interaction between the two components of the Fermi-Bose mix- 
ture is: 

/■ft/3 f 

S m M\^r^)=9BF J o dr J dr|V;(T,r)| 2 |0(T,r)| 2 , (4) 

where gs = 4^h 2 aB/mB and gsF = 2irh 2 aBF rnj = rriBm F / '{tub + rn F ), tub and m F 
are the masses of bosonic and fermionic atoms respectively, a# and as F are the s wave 
scattering lengths of boson-boson and boson-fermion interactions. 

Integral over Fermi fields is Gassian, we can calculate this integral and obtain the par- 
tition function of the Fermi system as a functional of Bose field 0(r, r). In the semiclassical 
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Thomas- Fermi approximation one has E] : 



pup 

Seff = / dTdvf eff {\4>{T,r)\), (5) 

JO 

o poo 

feff = -2 K ^J V~ede\n(l+e^) = 



00 e^de 



where e = p 2 /2m F , Jx = Hf — Vf(t) — ^ B jr|0(r, r)| 2 and « = 2 1/2 m F /2 /(37r 2 ft 3 ). 
So we can write the effective bosonic Hamiltonian in the form: 

H eff = j^|V0| 2 + (^(r)- MB )|0| 2 + 

+ fl0| 4 + /e//(|0|)}. (7) 

The first three terms in (jZj) have the conventional Gross-Pitaevskii form, and the last term 
is a result of boson-fermion interaction. In low temperature limit /2/(&bT) ^> 1 one can 
write f e ff(\(p\) in the form: 

/ e //(|0|) = -\^ 12 - ^ *(* B T) 2 /2 l/a - (8) 
As usual, //p can be determined from the equation 

N F = drn F (r), (9) 



where 

"fto = -K / - ■ .... . (10) 



2 J 1 + e^( £ -A) 

Let us consider now the 87 Rb and 40 K mixture with an attractive interaction between 
bosons and fermions 0. The parameters of the system are the following: a B = 5.25 nm, 
clbf = — 21.7i 4 'g nm. K and Rb atoms were prepared in the doubly polarized states 
\F = 9/2, mp = 9/2 > and |2, 2 >, respectively. The magnetic potential had an elongated 
symmetry, with harmonic oscillation frequencies for Rb atoms Ub,t = = 2tt x 215Hz 
and ujb, z = = 2tt x 16.3Hz,u>f = y m B j mF 1 ^ b ~ 1A7u>b 7 so that mB0J B /2 = 

mpojp/2 = Vq. The collapse was found for the following critical numbers of bosons and 
fermions: N Bc « 10 5 ; N K « 2 x 10 4 . 

At the zero temperature limit, expanding f e ff(\<j>\) up to the third order in g BB we 
obtain the effective Hamiltonian in the form: 

Heff = y"dr|^|V0| 2 + (y e //(r)-/i B )|0| 2 + 



+ !J f o 1 • -^g BF o 6 } . en; 



K 
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where 



Veff{r) = (l - ~«^ /J W) \rn B u 2 B (p 2 + AV) , (12) 

3 1/2 2 /in\ 

</ e // = 9b- tj W 3bf, (13) 

and p 2 = x 2 + y 2 . 

In principle, one can study the properties of a Bose-Fermi mixture with the help of 
feff © without any expansion. However, the form of the Hamiltonian gives the 
possibility to get a clear insight into the physics of the influence of the Fermi system on 
the Bose one (see discussion below). It may be easily verified that the expansion of the 
function f(x) = (1 + x) 5 ^ 2 (see Eq. (jHJ)) up to the third order in x gives a reasonably 
good approximation for f(x) even for rather large values of x, in contrast with the higher 
order expansions, so one can safely use Eq. (fTTj) as a starting point for the investigation 
of the properties of the Bose subsystem. 

In derivation of Eqs. (jllH13j) we also use the fact that due to the Pauli principle 
(quantum pressure) the radius of the Bose condensate is much less than the radius of the 
Fermi cloud Rp ~ a/ Jip/Vo, so one can use an expansions in powers of Vf(t)/ ' \ip- 

From Eq. (|T2*j) one can see that the interaction with Fermi gas leads to modification 
of the trapping potential. For the attractive fermion-boson interaction the system should 
behave as if it was confined in a magnetic trapping potential with larger frequencies than 
the actual ones, in agreement with experiment [H] . Boson-fermion interaction also induces 
the additional attraction between Bose atoms which does not depend on the sign oi g BF . 

The last term in H e ff (jllj) corresponds to the three-particle elastic collisions induced 
by the boson-fermion interaction. In contrast with inelastic 3-body collisions which result 
in the recombination and removing particles from the system [9 3 , this term for gsF < 
leads to increase of the gas density in the center of the trap in order to lower the total 
energy. The positive zero point energy and boson-boson repulsion energy (the first two 
terms in Eq. ()11J0 stabilize the system. However, if the central density grows too much, 
the kinetic energy and boson-boson repulsion are no longer able to prevent the collapse 
of the gas. Likewise the case of Bose condensate with attraction (see, for example, 0), 
the collapse is expected to occur when the number of particles in the condensate exceeds 
the critical value Nb c - 

The critical number N Bc can be calculated using the well-known ansatz for the Bosonic 
wave function: 

NbX V /2 / (P 2 + AV) \ 
w z aH m) exp \ 2w 2 a 2 )' 1 ' 

where w is a dimensionless variational parameter which fixes the width of the condensate 
and a = ^Jh/rriBUJB- 

In this case the variational energy Eb has the form: 

E B 2 + A1 ^ + cjNs + ^ 



N B hldB Aw 2 w 3 w 
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Figure 1: Variational energy jr—i — as a function of w for various numbers of bosons. 
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This energy is plotted in Fig.l as a function of w for several values of Nb- It is seen 
that when Nb < Nb c there is a local minimum of Eb which correspond to a metastable 
state of the system. This minimum arises due the competition between the positive first 
three terms in Eq. (|15j) and negative fourth term. The local minimum disappears when 
the number of bosons N B exceeds the critical value which can be calculated by requiring 
that the first and second derivatives of Eb vanish at the critical point. In this case the 
behavior of Eb is mainly determined by the second and fourth terms in Eq. (j!5|) . For 
Nk = 2 x 10 4 and clbf — ~ 19.44 nm we obtain Nb c ~ 9 x 10 4 in a good agreement 
with the experiment It is interesting to note that the critical number of Bose atoms in 
Bose-Fermi mixture is about two orders larger than the critical number for the condensate 
with a purely attractive interaction. For example, in the experiments with trapped 7 Li 
j2] it was found that the critical number of bosons is about 1000. 

Upon increasing the number of fermions, the repulsion between bosons decreases lead- 
ing to the collapse for the smaller numbers of the bosonic atoms. It should be also noted 
that an increase of the temperature results in a decrease of the local density of fermions 
and reduces the interaction energy between Bose and Fermi systems increasing the critical 
number of bosons. 

Finally, we make a short remark on the nature of the collapse transition. In this ar- 
ticle we found the instability point of the Bose-Fermi mixture with attractive interaction 
between components. A strong rise of density of bosons and fermions (in the collaps- 
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ing condensate enhances intrinsic inelastic processes, in particular, the recombination in 
3-body interatomic collisions, as is the case for the well-known 7 Li condensates [§]. How- 
ever, recently M. Yu. Kagan and coworkers suggested the new microscopic mechanism 
of removing atoms from the system which is specific for the Bose-Fermi mixtures with 
attraction between components and is based on the formation of the boson-fermion bound 
states [TU]. It seems that the description of the evolution of the collapsing condensate 
should include both these mechanisms. 
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